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We systematically study the BCS-BEC crossover and the quantum phase transition in ultracold 
^Li and ''"K atoms across a wide Feshbach resonance. The background scattering lengths for ®Li and 
''"K have opposite signs, which lead to very different behaviors for these two types of atoms. For ^"K, 
both the two-body and the many-body calculations show that the system always has two branches 
of solutions: one corresponds to a deeply bound molecule state; and the other, the one accessed by 
the current experiments, corresponds to a weakly bound state with population always dominantly 
in the open channel. For ®Li, there is only a unique solution with the standard crossover from the 
weakly bound Cooper pairs to the deeply bound molecules as one sweeps the magnetic field through 
the crossover region. Because of this difference, for the experimentally accessible state of '"^K, there 
is a quantum phase transition at zero temperature from the superfluid to the normal fermi gas at 
the positive detuning of the magnetic field where the s-wave scattering length passes its zero point. 
For ^Li, however, the system changes continuously across the zero point of the scattering length. 
For both types of atoms, we also give detailed comparison between the results from the two-channel 
and the single-channel model over the whole region of the magnetic field detuning. 



I. INTRODUCTION 

Ever since the experimental observation of the con- 
densation of fermion pairs in dilute gases of ultracold 
fcrmionic atoms, the study of the BCS-BEC crossover in 
ultracold fermionic gases has attracted much attention 
P, Q • The crossover in ultracold fermion gases is facil- 
itated experimentally by the Feshbach resonance, where 
the inter-atomic interaction can be tuned over a large 
range by varying the external magnetic field 0. Var- 
ious experiments have been done recently in ultracold 
fermionic gases in the vicinity of the Feshbach resonance 
to study the properties of the BCS-BEC crossover 

In most experiments, either an equal mixture of ^"^K 
atoms in the \F = 9/2, mp = -9/2) and the \F = 
9/2, mp = ^7/2) hyperfine states or an equal mixture 
of '^Li atoms in the \F = 1/2, mp = -1/2) and the \F = 
1/2, mp- = 1/2) states are prepared and magnetically 
tuned across a wide Feshbach resonance (Bq = 202G for 
-^^K and Bq = 834G for ^Li). Because of the nature of 
the wide resonance, for both atom species, their proper- 
ties close to the resonance point are very similar H, |^ |^ ■ 
However, they indeed have an important difference which 
is not quite emphasized before: the background scatter- 
ing lengths for and ^Li have opposite signs. This sign 
difference in the background scattering has profound ef- 
fects when the system is outside the unitary region. 

In this paper, we investigate the different properties of 
and ^Li atoms arising from their distinction in the 
background scattering. To describe the Feshbach reso- 
nance, we use the standard two-channel model which is 
characterized by three parameters determined from the 
scattering data: the resonance width^ the detuning, and 
the background scattering length [2 IE 0- The main 
differences between and ^Li are summarized as fol- 
lows: First, for the '^''K atoms, there always exist two 
stable solutions to the many-body gap and number equa- 



tions; while for the ^Li atoms, there is only one solution. 
The current experiments with ""^K atoms probe only one 
branch of the solution through the adiabatic sweep. For 
this experimentally relevant branch, the closed channel 
population is always limited to a small fraction, no mat- 
ter how far the magnetic field is detuned. This unusual 
property of ""^K has been noted before, in association 
with a somewhat different model Hamiltonian which is 
characterized by five experimental parameters |^ . For 
future experiments with '^''K, there might also be the 
possibility to switch to the other branch of the solution 
through a fast sweep, if the three-body collision does not 
turn out to be a significant obstacle to the stability of 
such a state (note that the three-body collision is not 
taken into account in the two-channel Hamiltonian). Sec- 
ond, a probably more interesting difference is that there 
exists a quantum phase transition at zero temperature for 
the '^''K atoms; while there is no such transition for ^Li. 
Note that here we limit ourselves to the case with equal 
populations for the two spin components, where there is 
no phase transition in the standard BCS-BEC crossover 
theory. Nevertheless, for the experimentally accessible 
branch of solution of the atoms, due to the effect 
of the background scattering in the two-channel model, 
a quantum phase transition from the superfluid state to 
the normal Fermi gas occurs on the BCS side of the res- 
onance, basically around the point where the scattering 
length goes across zero. For the ^Li atoms, however, 
when the scattering length crosses the zero point, which 
is on the BEC side, the change in the state of the system 
is completely continuous. 

In this work, we have also systematically compared the 
results from the two-channel model and from the single- 
channel model over the whole region of the magnetic field 
detuning (from the far negative to the far positive) . The 
standard single-channel model is characterized by only 
a single parameter, the scattering length at the corre- 
sponding field detuning. Our results extend some of the 
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comparisons in Refs. to the ofF-resonance region, with 
particular emphasis on the effects arising from the dif- 
ference in the background scattering of the two atom 
species. For easier comprehension of the properties of 
the many-body solutions, we also present the solutions 
for the two-body states with the same model Hamilto- 
nian for comparison. The paper is arranged as follows: In 
Sec. II, we present the model Hamiltonians, and discuss 
their two-body states for ^Li and ""^K atoms. Then, in 
Sec. Ill, we apply the two channel and the single channel 
models respectively to study the many-body properties of 
the ultracold atoms at zero temperature. The properties 
calculated include the energy gap, the chemical poten- 
tial, and the closed channel population. In Sec. IV, we 
analyze in detail the quantum phase transition for the 
^"K atoms on the BCS side. Finally, in Sec. V, we dis- 
cuss the finite temperature behavior, and calculate the 
critical temperatures and the pseudo-gaps corresponding 
to the different branches of solutions for ^Li and ^"K 
atoms. For this finite-temperature calculation, we follow 
the GqG method as reviewed in Ref. |^ to include the 
pair fluctuations. The results are summarized in Sec. VI. 



II. THE MODEL HAMILTONIANS AND THEIR 
TWO-BODY STATES 

The Feshbach resonance in ultracold atoms can be un- 
derstood on the grounds of two colliding fermions in dif- 
ferent hyperfine states. In dilute ultracold gases, the in- 
teraction can be described by the s-wave atomic scatter- 
ing. The magnetic field tunes the interaction strength by 
Zeeman shifting the energy level of the quasi-bound state 
(the Feshbach molecule state) relative to the continuum 
threshold of free atoms. The physics is well described by 
the following two-channel model 0, 



H 



cr,k q 

(U/V) ^ «k+q/2,T"-k-q/2,i"-k'-q/2,i»k'+q/2,T 



k,k',q 



g/^) K"-k+q/2a«k-q/2,T + ^.C.) , (1) 

k,q 



where aj^ ^ and fe^ are the creation operators for the 
open channel fermions (the atoms) and the closed channel 
bosons (the molecules), respectively; ek = Ti'k'^ /{2m) (m 
is the atom mass); V is the quantization volume, and a = 
t, i denotes the different hyperfine states of the atoms. 
This Hamiltonian is characterized by three parameters: 
the bare atom-molecule coupling rate the bare back- 
ground atom scattering rate J7, and the bare detuning 7. 
These three bare quantities are connected with the phys- 
ical ones gp,Up, through the standard renormalization 
relations: U = TUp , g — Tgp , 7 = 7p — Tg^/Uc, where 



^ and = - Ek 2^ 3- The physi- 



T^il + Up/U,j ..... =-^k% 
cal quantities gp, Up,'jp are determined from the scatter- 



ing data as Up = Airf^'abg/m^ g^ = Anh^abgAB^co/iTT-, 
and 7p = fico{B — Bq) {^co is the difference in magnetic 
moments between the two channels) jlOj . where we have 
assumed that the s-wave scattering length near resonance 

has the form as — a^g ~ Bq ) ' ^^^^ ^^a back- 
ground scattering length, AB as the resonance width, 
and i?o as the resonance point. 

If the population in the closed channel is negligible, one 
can adiabatically eliminate the molecular modes 6q in the 
Hamiltonian (1), and arrive at the following simplified 
single-channel model 



(2) 
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+ {V^f /V) Y «I+q/2.T«-k-q/2,i"-^k'-q/2,i«k'+q/2,T, 
k,k',q 

where V^^^ is connected with the physical one 
Upf^ by the renormalization relation V^^^ = 

U^f^l + U^ff/Uc)~\ and U^ff is determined solely 
by the scattering length as U^^^ — Airh^as/ra 
at the corresponding magnetic field detuning. So, a 
single-channel model does not explicitly depend on the 
resonance width or the background scattering rate, 
it only depends on their combination through the 
expression of a^. In this sense, a single- channel model 
is simpler and more universal. However, given a new 
configuration, we typically do not know whether the 
closed channel population is negligible prior to the 
calculation, therefore the validity of the single-channel 
model is not self-evident. We will have some systematic 
comparison of the results from the two-channel and the 
single-channel calculations for and ^Li over the 

whole region of the magnetic field detuning. Similar 
comparison has been done before in the unitary region 
0,1a], where the background scattering can be neglected. 

To have some intuitive understandings of the Hamil- 
tonian above, we first look at its bound states for two 
particles. For a two channel model, the eigenstate of the 
Hamiltonian in the case of two particles can be written as 
the superposition of the molecular bosons in the closed 
channel and the fermion pairs in the open channel: 



I*) 



ka 



Cbbl |0) , 



(3) 



where Ck and Ch are the superposition coefficients, and we 
have assumed that the center-of-mass momentum q = 
as it is decoupled from the relative momentum in the two- 
body case. The Schrodinger's equation = E\^>) 
then gives: 
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where we have used the renormahzation relation 

1/ [u, - 9lh,] = 1/ [c/ - gVi] + (i/v) Ek [1/ (2^k)]. 

We solve the bound state energy levels of ^Li and 
^"K, respectively. The parameters for the ^Li atoms are 
abg = — 1405ao, AJ5 = — 300G and /ico ^ '^^J-b, where 
ao is the Bohr radius and /is is the Bohr magnetic mo- 
ment. The parameters for the ^"^K atoms are a^g = 174ao, 
ABq = 7.8G, and fico ~ 1.68/i_B- The results of our cal- 
culation are plotted in Fig. 1. The outstanding feature of 
the figures is that there are two branches of bound state 
for "^''K (Fig. 1(a)). One branch is a deep-bound molecule 
state which never crosses the continuum threshold (zero 
energy in the plots). In this state, the closed channel 
fraction (|cb| ) becomes dominant as one approaches the 
negative detuning. The other branch (the upper one) is 
a weakly bound pair state which crosses the continuum 
threshold at 7 = 0. In this state, the closed channel 
fraction is always limited to a small fraction (< 6%), 
no matter how much the magnetic field is detuned. In 
experiments, one approaches the bound state by adiabat- 
ically tuning the magnetic field, starting from the atom 
continuum. Therefore the branch of state with weakly 
bound pairs is the one accessed in the current experi- 
ments. If one wants to go to the deeply bound molecule 
state, one either needs a fast magnetic field sweep with 
a rate larger then the energy splitting between the two 
bound states, or use some rf pulses to couple these two 
branches. For the ^Li atoms, in contrast, there is only 
one branch of solution, which merges into the continuum 
at 7 = and approaches the deep-bound molecule state 
as one increases the negative field detuning. This differ- 
ence between the two atom species comes from the differ- 
ent natures of their background scattering interactions. 
For ''°K with a positive background scattering length, 
as has been discussed in Refs. |3, |9|, (the Hamilto- 
nians therein are characterized by five parameters and 
are somewhat different from the model here), there is a 
weakly bound state in the open collision channel. The 
avoided level crossing between this weakly bound state 
and the Fcshbach molecule level gives the two branches 
of solutions. For ^Li with a negative background scatter- 
ing length, there is no bound state in the open channel, 
and thus there is only a single branch of solution. 

For comparison, we have also calculated the two- 
body bound states using the single channel Hamiltonian 
(2). In the single-channel approach, one simply replaces 

Up - in Eq. (4) with U^^f, so the left side of 

Eq. (4) becomes energy independent. Equation (4) then 
gives the well known expression of the bound energy 
E = —h^/ {ma^) in the region Os > for the single- 
channel mode. The results of the calculation are shown 
as the dash-dotted curves in Fig. 1. The single channel 
calculation gives almost the same bound state energies 
for the weakly bound state of atoms (it neglects the 
other branch of solution); while it deviates from the re- 
sults of the two channel calculation as the atomic scat- 
tering length as approaches its zero point (corresponding 
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FIG. 1; The bound state energy levels for two colliding 
atoms of (a)'*''K, (b)''Li. The solid curves are the deeply 
bound states, the dashed curve is the weakly bound state 
in ''"K, and the dash-dotted curves are the results from 
single channel Hamiltonian. To compare with the many- 
body calculations that follow, the unit of energy is taken 
to be Ef = ?i^(37r^n)^/^/2m, where the total particle den- 
sities (n) are chosen to be typical values in the experiments 
(n ~ 1.8 X lO^^cm"^ for '^°K, and n ~ 2.9 x 10^^ cm~^ for ''Li 
JJ )• 

toB-Bo^ 7.8Gfor "Ok, and B - Bq = -300G for ^Li). 

III. THE MANY-BODY STATES 

To find out the many-body state at zero temperature, 
one needs to minimize the energy corresponding to the 
Hamiltonian (1) under the constraint of total particle 
number conservation N = X^o- k ''k o-'^k.o- + ^^^fo^foq. 
One therefore minimizes = {H — /iiV) , where /i is the 
Lagrange multiplier with the physical meaning of the 
chemical potential. We follow the standard BCS-BEC 
crossover theory j which assumes a mean field for the 
molecule and the pair operators, with (6q) = (6o)(5qo = 

~ J2u ('^-k.i'^kj) / (7 — 2/Li), where the second 

equality comes from the Heisenberg equation for the op- 
erator 6q. Under such an assumption, one can find out 
the explicit expression of the energy f2, whose extrema 
conditions — ^ — ^ yield the standard gap equa- 

tion and the number equation: 

1 _ 1 ^ 1 1 

ntot = 2"^. + ^E(l-^^)' (6) 
^ k 
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where the quasi-particle excitation energy E]^ = 

^(ek-/^)2 + |Ap, and the gap A ^ g{bo)/VV + 

iU/V)J2k{a-u,iaKl) = zi,{bQ)/VV {zb = g ~ 
U il ~ 2n) /g ^ gp - Up {jp - 2/i) /gp). In Eq. (6), ritot 
denotes the total particle density, and — |(&o)l^ /"l^ is 
the density of the closed channel molecules. 

With Eqs.(5) (6), we can solve for the gap, the chemi- 
cal potential and the closed channel fraction at different 
magnetic field detunings. The results of the calculations 
for ""^K and for ^Li are plotted in Fig. 2 and Fig. 3, 
respectively. For ^°K, similar to the two-body case, Eqs. 
(5) (6) support two branches of solutions. The proper- 
ties of the first branch of solutions is shown in Fig. 2 
(a-c). Compared with the results in the previous sec- 
tion, it is easy to see that this solution corresponds to 
the deep-bound state in the two-body energy spectrum. 
The closed channel fraction increases from in the BCS 
limit to 1 in the BEC limit, and the maximum in the 
gap appears on the BCS side of the resonance, near the 
Us = point. As the chemical potential always remains 
negative for this branch, the system is essentially BEC- 
like at all magnetic field detunings. The properties of 
the second branch of solution is shown in Fig. 2 (d- 
f), where the closed channel fraction is non-monotonic, 
peaks around B — Bq ^ —5.2G with a fraction ~ 5.1%, 
in rough agreement with the calculations in Refs. 
where a different model is used. The chemical poten- 
tial in this branch crosses zero on the BEC side of the 
resonance, demonstrating the existence of a BCS-BEC 
crossover for this branch. The most intriguing feature of 
the second branch is that there is a quantum phase tran- 
sition on the BCS side of the resonance, beyond which the 
superfluid disappears, and the system phase-transits into 
a normal fermi gas. We will discuss this phase transition 
in more details in the next section. The properties of ^Li 
atoms are shown in Fig. 3. The closed channel fraction 
remains small throughout the crossover region, and only 
becomes appreciable deep in the BEC region. The peak 
in the gap occurs again near the as — point (now on the 
BEC side). The chemical potential crosses zero around 
B- Bo 85G, denoting a BCS-BEC crossover. 

For comparison, we have calculated the gap and the 
chemical potential using the single channel Hamilto- 
nian. For the single-channel approach, one just replaces 

Up ~ ^ ^2fi "^ith Upf^ in the gap equation (5) and drops 
the closed channel population contribution 2nh in the 
number equation (6). The results are shown in Fig. 2 
and Fig. 3 in dash-dotted curves. The results from the 
single channel calculation should coincide with the two 
channel calculation provided that the closed channel pop- 
ulation is small. From Fig. 2 and Fig. 3, it is obvious 
that, except for the deep-bound state solution of '^''K, 
the single channel Hamiltonian reproduces the results of 
the two channel Hamiltonian for ^"K and ^Li for most of 
the crossover region. Notably, the single channel descrip- 
tion is a fairly good approximation for almost all mag- 
netic field detunings for the weakly bound state solution 
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FIG. 2: Results from the manybody calculation for '"^K atoms 
at zero temperature, (a-c) The gap, the closed channel frac- 
tion, and the chemical potential of the solution which corre- 
sponds to the deeply bound state in the twobody case; (d-f) 
The gap, the closed channel fraction, and the chemical poten- 
tial of the solution which corresponds to the weakly bound 
state in the twobody case. The solid curves are the results 
using two-channel Hamiltonian. The dash-dotted curves are 
the results using single channel Hamiltonian. The unit of en- 
ergy is the same as defined in the caption of Fig. 1. 



of ^''K, as the closed channel fraction in this branch re- 
mains small. In both cases, the single channel description 
breaks down near the point where as — 0. 



IV. QUANTUM PHASE TRANSITION FOR 
ATOMS 



To see more clearly the quantum phase transition, we 
plot the gap and the closed channel fraction from the two 
channel calculation near the phase transition point on a 
logarithmic scale in Fig. 4. Both the gap and the closed 
channel fraction show a precipitous drop at the magnetic 
detuning B~Bq ~ 7.8G, the point where — 0. Beyond 
that point, the gap equation and the number equation do 
not have solutions (for the upper branch). This is a clear 
signature of a quantum phase transition. The ground 
state of the '^^K atoms changes from a BCS superfluid 
state into a normal state with the gap A = through a 
second order phase transition. In contrast, for ^Li atoms 
(Fig. 3), there is just a smooth crossover for both the 
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FIG. 3: Results from the manybody calculation for "Li atoms 
at zero temperature. (a)The gap and the closed channel frac- 
tion (inset); (b) The chemical potential. The solid curves are 
the results using two-channel Hamiltonian. The dash-dotted 
curves are the results using single channel Hamiltonian. 



gap and the closed channel population at the point where 
a,^0{B-Bo^ -300G). 

The origin of the phase transition comes from the re- 
pulsive background inter-atomic interaction between the 
^•^K atoms in the two different hyperfine states. The re- 
pulsive interaction does not support a superfluid state 
which can be adiabatically connected to the BCS super- 
fluid in the region < (although the system indeed has 
a superfluid solution corresponding to the deep-bound 
state, it is not adiabatically connected). So, across the 
point Os = 0, there must be a phase transition, and the 
atoms transit into a normal phase when > 0. If one 
starts from this normal Fermi g the magnetic field 
decreases towards the resonance {Bq = 202G), the atoms 
will first undergo a second order phase transition and be- 
come a BCS superfluid at B - Bq 7.8G. If the field 
continues to decrease, the chemical potential will even- 
tually become negative at S — So ^ — 0.34G, and the 
atoms enter the EEC region with a smooth crossover. 



V. THE CRITICAL TEMPERATURES FOR 
DIFFERENT BOUND STATES 

To characterize the crossover properties of different 
bound states for different atom species, we calculate the 
critical temperature at which all the pairs become non- 
condensed. At finite temperature T, instead of the energy 
minimization, one needs to minimize the thermodynam- 
ical potential = — Tln[tr(e~^/"^)]. If one takes the 
mean- field approach outlined in Sec. II, from the ex- 



FIG. 4: The enlarged (a) gap and (b) closed channel fraction 
as a function of magnetic field detuning for atoms near 
the phase transition {B ~ Bo ^ 7.8G) (corresponding to the 
weakly bound branch of solution). 



trema conditions 



dn _ dn 



— ^ — 0, one gets the following 



9 (bo) ^^l 

finite-temperature version of the gap equation and the 
number equation |^ll3| 



1 



9i 



2Eu 



2ek 



ritot = 2nf, 



Eu 



where f{x) = l/(e^/'"'^^ + 1) is the fermi distribution, 
and ks is the Boltzmann constant. 

This simple mean- field approach however, is not ade- 
quate for the calculation of the critical temperature ex- 
cept in the case of the BCS limit. The reason is that at 
finite T, due to the thermal fluctuations, there exist non- 
condensed pairs (superposition of fermion pairs and non- 
condensed bosons), in addition to the condensed fermion 
pairs and the condensed molecular bosons Hia. The 
non-condensed pairs also contribute to the gap for the 
fermionic excitations. Accommodating this contribution, 
the gap is replaced by the total gap |A|^ = |As|^-f| Apg \ ^ , 
where As — Zb{bo)/VV is the same as the gap defined 
in the previous section (the superfluid order parameter), 
and the pseudo-gap Apg comes from the contributions of 
the non-condensed pairs 0,0]. With this simple replace- 
ment, the above gap equation and the number equation 
from the mean-field approach remain valid. One just 
needs to interpret A as the total gap, and thus the emer- 
gence of A does not correspond to the critical transition 
between the normal and the superfluid phase. Instead, 
the gap shows up at a temperature T*, below which the 
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non-condensed pairs are formed. As one further decreases 
the temperature, these pre-formed pairs finally condense 
and become a superfluid at the critical temperature Tc 
with Tc < T*. The gap at the critical temperature Tc 
solely comes from the pseudo-gap Apg, contributed by 
the non-condensed pairs (as = at that point). To 
calculate this phase transition temperature Tc, one needs 
to breakup the total gap A into the superfluid order 
parameter Ag and the pseudo-gap Apg, which requires 
knowledge of the dispersion relation of the pair excita- 
tions. As shown in Ref. 01, the dispersion relation for the 
pair excitations is still quadratic with a self-consistently 
determined effective mass M*. With that formalism, we 
calculate the critical temperature Tc and the pseudo-gap 
at Tc for different branches of states in and ^Li over 
the whole region of the magnetic field detuning. 

The results of the calculation are plotted in Fig. 5 
for atoms, and in Fig. 6 for ^Li atoms. Fig. 5 (a) 

(b) correspond to the deeply bound state, and Fig. 5 

(c) (d) correspond to the case with weakly bound state. 
It is clear from the Tc calculation that the BCS-BEC 
crossover picture is evident only for the second branch of 
the solution (Fig. 5 (c) (d)). For this branch, the criti- 
cal temperature becomes exponentially small on the BCS 
side and almost constant on the BEC side as one expects. 
For the other branch of the solution corresponding to the 
deep-bound state, the critical temperature only varies in 
a very small range, and converges to a constant value on 
either side of the resonance. One can also see that except 
in the BCS limit, the pseudo-gap is still significant at the 
critical temperature. 
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FIG. 5: The critical temperature Tc and the pseudogap at 
the critical temperature Ac as a function of magnetic field 
detuning for the two solutions in ^"K atoms ((a)(b) corre- 
spond to the deeply bound branch, while (c) (d) correspond 
to the weakly bound branch). Tc is in the unit of Tf, the 
Fermi temperature associated with the Fermi energy defined 
in the caption of Fig. 1. 



VI. SUMMARY 

We have systematically studied the BCS-BEC 
crossover and the quantum phase transition in ultracold 
^Li and ""^K atoms across a wide Feshbach resonance. As 
the background scattering lengths have opposite signs for 
these two types of atoms, their properties are quite differ- 
ent outside of the near resonance region. Both the two- 
body and the many-body calculations show that for "'"K, 
there always exist two branches of solutions, correspond- 
ing respectively to the deep-bound state and the weakly 
bound states. The latter one with weakly bound pairs is 
the branch that is accessed by the current experiments 
through adiabatic sweep. For this branch of solution, 
there is a quantum phase transition on the BCS side of 
the resonance, where the atoms phase-transit from the 
superfiuid state to the normal fermi gas as the scattering 
length crosses zero. This kind of phase transition is ab- 
sent in the ^Li atoms, where the change in the state of 
the system is completely continuous when the scattering 
length goes across its zero point. 

The calculations here emphasize the importance of 
the background scattering length outside of the near- 
resonance region, and the different properties of the ^Li 
and "^^K atoms. In this sense, ^Li and ''"K are good rep- 
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FIG. 6: The critical temperature Tc and the pseudogap at 
the critical temperature Ac as a function of magnetic field 
detuning in ^Li atoms. Tc is in the unit of Tp, the Fermi 
temperature associated with the Fermi energy defined in the 
caption of Fig. 1. 
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resentatives of two kinds of atoms which have different 
background scattering properties. The calculations here 
may also have some interesting experimental indications. 
For instance, one can verify the quantum phase transi- 
tion for the atoms on the BCS side of the resonance, 
and the continuous crossover for the ^Li atoms on the 
BEC side of the resonance, although in both cases the 
scattering length crosses its zero point in a similar way. 



For the K atoms, there might also be the possibility 
to access the deeply bound branch of solution in future 
experiments through either a fast magnetic field sweep 
or some rf transitions. 

This work was supported by the NSF awards 
(0431476), the ARDA under ARO contracts, and the A. 
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discussions. 
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